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ABSTRACT 

We performed a detailed analysis of the main theoretical uncertainties affecting the age 
at the lithium depletion boundary (LDB). To do that we computed almost 12 000 pre- 
main sequence models with mass in the range [0.06, 0.4] Mq by varying input physics 
(nuclear reaction cross-sections, plasma electron screening, outer boundary conditions, 
equation of state, and radiative opacity), initial chemical elements abundances (total 
metallicity, helium and deuterium abundances, and heavy elements mixture), and 
convection efficiency (mixing length parameter, aMu)- As a first step, we studied the 
effect of varying these quantities individually within their extreme values. Then, we 
analysed the impact of simultaneously perturbing the main input/parameters without 
an a priori assumption of independence. Such an approach allowed us to build for 
the first time the cumulative error stripe, which defines the edges of the maximum 
uncertainty region in the theoretical LDB age. 

We found that the cumulative error stripe is asymmetric and dependent on the 
adopted mixing length value. For aML = 1-00, the positive relative age error ranges 
from 5 to 15%, while for solar-calibrated mixing length, the uncertainty reduces to 
5 — 10%. A large fraction of such an error (« 40%) is due to the uncertainty in the 
adopted initial chemical elements abundances. 

Key words: Methods: numerical - stars: abundances - stars: evolution - stars: 
fundamental parameters - stars: low-mass - stars: pre-main sequence 


1 INTRODUCTION 

Lithium burning via the proton capture ^Li(p, a)a becomes 
efficient in stellar conditions as temperature rises above 
2.5 X 10® K. In the mass range 0.06-0.4 Mq (the exact val¬ 
ues depending on the chemical composition), such a tem¬ 
perature is already reached during the pre-main sequence 
(pre-MS) phase allowing lithium to be completely destroyed 
in fully convective structures. Since the larger the mass and 
the higher the rate of central temperature increase and, con¬ 
sequently, the earlier the onset of lithium burning, the age 
at which lithium is fully depleted strongly depends on mass. 
Thus, in stellar associations or clusters with ages between 15 
and 350 Myr, one would expect to observe a sharp transition 
in the very low mass regime between stars with and without 
photospheric lithium. Such a transition, usually called lith- 
ium depletion boundary (LD B), is an age indicator (see e.g. 
iD’Antona fc Mazzitellilll994l L 
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The age at the LDB can be provided only by theoret¬ 
ical stellar evolutionary models. The accuracy of their pre¬ 
dictions relies on the adopted input physics (e.g. equation 
of state, radiative opacity, nuclear reaction cross-sections, 
etc.), chemical composition (initial metallicity, helium and 
deuterium abundance, heavy element mixture, etc.), and the 
numerical scheme describing macroscopic physical processes 
(convection, rotation, etc...). All these ingredients are af¬ 
fected by not negligible uncertainties, which translate into a 
theoretical error in the LDB age estimate. 

Given the primary importance of evaluating the age of 
young stellar clusters, several studies have been devoted to 
quantify the main theoretical and observational uncertain¬ 
ties affecting the LDB method. A first idea of the theoretical 
uncertainty is provided b y comparing models co mputed by 
different groups (see e.g. I Jeffries fc Navlonl200lll . However 
such an approach is neither able to quantify the individual 
uncertainty sources nor able to fully exploit the whole the¬ 
oretical uncertainty, as several input physics are in common 
among the different evolutionary codes. 

A better and more detailed way to proceed is to vary 
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a single ingredient at a time and to study its impact on 
the LDB age (see e.g. Bildsten_e£^ 1997|, lUshomirskv et ^ 
1 19981 . iBurke. Pinsonneault fc Sillsl 2004). The weakness of 


this approach is that it is unable to quantify the interactions 
between the different uncertainty sources. 

A more robust procedure, although much more compu¬ 
tationally expensive, consists in varying simultaneously the 
input physics, parameters, and chemical abundances adop¬ 
ted in stellar computations within their current uncertainty 
range. Such an estimate of the theoretical uncertainty af¬ 
fecting the LDB ages is still lacking and, thus, we decided 
to provide it. We addressed also for the hrst time the ana¬ 
lysis of the impact of the current uncertainty in the adop¬ 
ted chemical composition on the LDB age estimates. To do 
that, we followed the same procedure adopted in our previ¬ 
ous studies on the cumulative physical uncertainty affecting 
low-mass star models from the m ain sequence to the He- 
burning phase (IValle et al.l 12013^1^ . that is, a systematic 
and simultaneous variation of the main input physics on a 
fixed grid. 

The number of models required to fully cover the whole 
parameter space is huge. For this reason, as a first step, we 
started to study the effect of individual input perturbations 
to discriminate between negligible and not negligible error 
sources. In Sect. [3] we computed sets of perturbed models 
by changing one input physics at a time keeping fixed the 
chemical composition and in Sect.|T]we computed sets of per¬ 
turbed models by varying individually the initial metallicity, 
helium and deuterium abundances, and heavy elements mix¬ 
ture at a time at fixed input physics. Finally, in Sect. [S] we 
computed sets of models by simultaneously varying all the 
ingredients that have a not negligible effect on the LDB age. 
Section reports the main conclusions. This work required 
the computation of about 12 000 pre-MS tracks. 


2 THE REFERENCE SET OF MODELS 

We computed the models with the most recent release of the 
PROSECCO (Pisa Raphson-NewtOn Stellar Evolution 
Computation COde) stellar e volution code derived from 
the FRANEC code (see e.g. iDegl’Innocenti et al.l l2008ll . 
A detail ed description of the standard models can b e 
found in Tognelli, Pra da Moroni fc Degl’Innocentil ll201llf . 
Toen elli. Degl’I nnoce nti fc Prada Mororii ll2012lf and 
DeU’Qmodarme et al.n2012h : here we limit the discussion 
only to the variations with respect to the previous version. 

The reference models are standard pre-MS tracks, which 
do not take into account accretion, rotation, and magnetic 
fields. The evolution starts from a fully convective and ho¬ 
mogeneous model on the Hayashi track, with a central tem¬ 
perature Tc ~ 10® K, low enough for deuterium burning to 
be completely inefficient. Although the initial model is not 
physically realistic, as it does not result from a self consist¬ 
ent protostar evolution, we checked that varying the initial 
central temperature (hence radius) in the range [10®, 10®] K 
does not affect the LDB. 

We adopted the 2 006 release of the OPAL equation of 
state (OPAL06 EOS; [Rogers fc Navfonovll2003 l. However, 
in order to compute objects less massive than 0.1 Mq, in 
the current version of the code we included also the EOS by 


Saumon. Chabrier fc van Horn! ll 19951 . SCVH95), for tem¬ 
peratures and densities not covered by the OPAL06 EOS. 

We used t he sa me radiative (OPAL 2 005 and 
^rgusonjt_^y__ 200^)__and co nduct ive (jPotekhini 
19991. [Sh ternin fc Yako vlev! l2006lj o paci ties as in 
Tognelli. Prada Moroni fc Degl’Innocentil 1 I 2 OIIII but 


for a different solar heavy elements mixture, namely the 
lAsplund et al.l (l2009l . AS09). 

Nuclear reaction rates relevant for the pre-MS evol- 
ution have been taken from the NACRE compilation 
(Angulo^j^^ 19^), with the exception of p(p,e^z/)^H 


( Marcucci. Schiavilla fc Vivianill2013l . 


2H(p,7)®He 
and ^H(^H , n)®He 


(Descouvemont et al 


Tognelli et al.l[2015l L 

2004h. ^H(%.pi®H 


( Tumino et al.l l2[)14^.®Li('p,®He) Q 
llLamia et al.l l2013lj . and ^Li(p,a)a (iLamia et al.l l20ld L 
Nuclear reaction rates between bare nuclei have been correc¬ 
ted to account for plasma elec tron screening. We implemen- 
ted the weak ([S alpeteil |l954| ), weak—intermediate—stron g 
dPewitt. Graboske fc Cooperl 19731. iGraboske et al. 


and strong ( Itoh. Totsuii &: Ichimarrilll977l . Iltoh et al.lll9'^ 

screening. Then, an adaptive procedure determines the 
most suitable screening factor to be used in each mass-shell 
of the model. For what concerns the burning of the light 
elements relevant for this work (mainly and ^Li), the 
intermediate/strong screening is generally adopted. 

Outer boundary conditions (BCs) have been 
obtained by the detaile d atmospheric models by 
Brott fc Hauschildd ll 2005l. BH05) a s discussed in 

Tognelli. Prada Moroni fc Degl’Innocentil ll201llf . Con- 


vection is treated ac cording to the mixing length theory 
([BohnnVitensd [195^ following the formalism described 


Cox fc Giulil i 19681 1. Usually the solar-calibrated mix¬ 


ing length value (oml) is adopted for different stellar 
masses and evolutionary phases. However, there are 
no compel ling theoretical arguments for s uch a choice 
(see e.g., iLudwig. Frevtag fc Steffeiil 1 19991 . iTrampedachl 
l2007tl . In particular, there are several hints that in 
pre-MS phase the supera diabatic convection might be 
less efficient (see e.g. iD’Antona &: MontalbanI 2003, 


Stassunet^ 2004. Mathieu et al.l 2007 . StassunI 


20081 . 


Somers fc Pinsonneault! l2014l. and references therein). 


Using our models a good agreement with eclipsing binaries 
and surface lithium ab undances in young clusters is ob¬ 
tained with aML^ 1 ( ^nn aro. Prada Moroni fc^Tognellil 
I2OI2I . [Tognelli. Degl’Innocenti fc Prada Moronil I2OI2I I. We 
adopted aML= 1-00 in our reference models. For the sake of 
completeness, in the following we will show also the effect 
of using our solar-calibrated value aML= 1- 743 

The reference set of models have been computed ad¬ 
opting the solar heavy element mixture by lAsplund et al.l 
( 2009!) and a heliu m-to-metal enrichment ratio AT/AZ = 2 
I Casagrandell200^ . With this choices, [Fe/H] = -1-0.0 trans¬ 
lates into an initial helium abundance Y = 0.274 and 
a global metallicity Z = 0.013 (see Sect. 11. The initial 


^ The solar model has been computed using an iterative pro¬ 
cedure to adjust the initial helium abundance, metallicity, and 
mixing length parameter in our 1 Mq stellar model, in order to 
reproduce, within a given numerical tolerance (< 10““^), at the 
age of the Sun (4.57 Gyr) its observed radius, luminosity, and 
(2/A)ph. 
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Figure 1. Age at which lithium is depleted by a factor of 100 
(LDB) as a function of luminosity for our reference set of models 
with Z = 0.0130, Y = 0.274, Xa = 2 X 10“®, and aya.= l-OO. 
Stellar masses (in Mq) are overplotted. 

abundances of the light elements are the sam e used in 
iTognelli. Degl’Innocenti fc Prada Moronil ll2012l ~). 

2.1 Lithium depletion boundary 

We computed 36 pre-MS tracks in the mass range [0.06, 0.40] 
Mq , in which lithium is completely destroyed in a fully con¬ 
vective object. We focused on this mass range because for 
the chosen value [Fe/H] = +0.0 objects less massive than 
about 0.06 M© never reach the temperature required to des¬ 
troy lithium, while for M > 0.4 Mq a radiative core develops 
before lithium is completely depleted. 

We defined the LDB as the model at which lithium 
abundance in mass is reduced by a factor of 100 with 
respect to the original one. In order to compute a re¬ 
liable lithium depletion evolution and, consequently LDB 
ages, the time-step adopted in model computations must 
be chosen with care. A time-step simply tuned to obtain 
an accurate pre-MS evolution would result in a too crude 
time resolution and in an underestimate of the lithium de¬ 
pletion at a given age, which in turn would lead to an 
overestimate of the LDB age dPiau &: Turck-Chiezd l2002l . 

iBurke, Pinsonneault fc Silla 2004l l. This is the consequence 

of ^Li-burning occurring over temporal-scales much shorter 
than the pre-MS evolutionary ones. Being aware of this, 
in the routine that provides the time-step for each itera¬ 
tion we introduced an additional condition that constraints 
the time-step to be short enough to lead to a relative lith¬ 
ium abundance variation lower than 0.1%. Moreover, such a 
choice allows us to obtain also a good numerical resolution 
in the ^Li abundance to precisely identify the LDB point. 

Fig.m shows our reference LDB curve in the (logL/L©, 
age) plane obtained from the reference set of models. The 
complete set is listed in Table [1] with the corresponding age, 
log L/L©, logLeff, radius, and surface gravity values at the 
LDB for each masfl 

We also provide a fit of the log age[yr] - log L/L© curve 


Table 1. LDB quantities for the reference set of models with 
Z = 0.0130, Y = 0.274, Xa = 2 x 10“®, and oml= 1-00. 


Aftot 

(M©) 

Age 

(Myr) 

log L/L© 

log i//R© 

logg 

(cm s“^) 

log Teff 
(K) 

0.060 

359 

-3.555 

-0.972 

5.160 

3.359 

0.061 

300 

-3.449 

-0.954 

5.132 

3.377 

0.062 

260 

-3.364 

-0.939 

5.108 

3.390 

0.063 

232 

-3.298 

-0.925 

5.088 

3.400 

0.064 

211 

-3.243 

-0.913 

5.070 

3.407 

0.065 

194 

-3.195 

-0.902 

5.055 

3.414 

0.066 

180 

-3.149 

-0.891 

5.040 

3.420 

0.068 

158 

-3.069 

-0.872 

5.015 

3.430 

0.070 

142 

-3.002 

-0.855 

4.993 

3.439 

0.072 

129 

-2.941 

-0.839 

4.972 

3.446 

0.074 

119 

-2.888 

-0.824 

4.954 

3.452 

0.075 

115 

-2.864 

-0.817 

4.946 

3.454 

0.076 

111 

-2.840 

-0.809 

4.937 

3.456 

0.077 

107 

-2.817 

-0.803 

4.929 

3.459 

0.078 

104 

-2.797 

-0.796 

4.922 

3.461 

0.079 

101 

-2.776 

-0.790 

4.915 

3.462 

0.080 

98 

-2.757 

-0.783 

4.908 

3.464 

0.082 

93 

-2.720 

-0.771 

4.895 

3.467 

0.084 

88 

-2.686 

-0.760 

4.882 

3.470 

0.086 

84 

-2.653 

-0.749 

4.869 

3.473 

0.088 

81 

-2.623 

-0.738 

4.858 

3.475 

0.090 

78 

-2.595 

-0.728 

4.847 

3.477 

0.095 

71 

-2.529 

-0.703 

4.822 

3.481 

0.100 

65 

-2.469 

-0.681 

4.799 

3.485 

0.110 

57 

-2.367 

-0.640 

4.759 

3.490 

0.120 

50 

-2.277 

-0.603 

4.723 

3.494 

0.130 

45 

-2.195 

-0.570 

4.691 

3.498 

0.140 

41 

-2.120 

-0.539 

4.662 

3.501 

0.160 

35 

-1.984 

-0.484 

4.611 

3.508 

0.180 

31 

-1.867 

-0.437 

4.568 

3.514 

0.200 

28 

-1.767 

-0.397 

4.532 

3.518 

0.225 

24 

-1.657 

-0.352 

4.493 

3.523 

0.250 

22 

-1.560 

-0.312 

4.459 

3.527 

0.300 

18 

-1.392 

-0.243 

4.401 

3.535 

0.350 

15 

-1.248 

-0.185 

4.353 

3.542 

0.400 

13 

-1.124 

-0.136 

4.312 

3.549 


Table 2. Coefficients of the third-order polynomial fit of the LDB 
age versus luminosity. 


ao 

ai 

a2 

as 

6.6408 

-0.46379 

-0.031392 

-0.014906 


for the reference set of models using a third-order polynomial 

logage[yr] = ao + ai x log L/L© + a2 x (logL/L©)^ 
+03 X (log L/L©)® 

The fit parameters oq, fli, 02, and 03 are given in Table [2] 
the mean accuracy of the fit is better than 1% over the whole 
selected luminosity range. 


^ Table[T]along with LDB ages for several [Fe/H] values is available 
in electronic form at the url: 

http://astro.df.unipi.it/stellar-models/ldb 


‘log’ stands for the base-10 logarithm. 
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3 INDIVIDUAL INPUT PHYSICS 

UNCERTAINTIES 

The results of stellar model computations depend on the ad¬ 
opted input physics, such as nuclear reaction cross-sections, 
plasma electron screening, radiative opacity, equation of 
state (EOS), and DCs. All these ingredients are still affected 
by a not negligible uncertainty which directl y translates 
into a n uncertainty in stellar model outcomes ijValle et alJ 
l2013all^ . In this section, we focus only on the input physics 
that might affect the lithium burning in pre-MS of very low 
mass stars. 

As a first estimate of the physical uncertainty impact 
on the LDB age, we varied a si ngle input physics at a time 
keepi ng all the others fixed (see lBurke, Pinsonneault fc Silld 
|2004 for a similar investigation). More in detail, when the 
error on a given input physics was available, we computed 
two additional sets of perturbed models by adopting respect¬ 
ively the highest and lowest value of the input physics given 
by its uncertainty. Then, the LDB ages provided by the per¬ 
turbed models have been compared with the reference one. 
Unfortunately, not for all the analysed input physics an un¬ 
certainty evaluation is available; this is for example the case 
of the EOS and outer BCs. Each of these cases has been con¬ 
veniently treated, as discussed in the following sub-sections, 
by substituting the adopted tables. Table|3]lists the analysed 
input physics with the related assumed uncertainty/range 
of variation, when present, or the alternative input physics. 
Note that, only in the case of the nuclear reaction cross- 
sections the listed errors are la, while in the other cases 
they represent the extreme values of the variability region. 

Where not explicitly stated, all the models have been 
computed for the reference chemical composition and mixing 
length parameter, as described in Sect. [5] 

Some of the cases analysed in the following sub-sections 
have b een already discussed in iBurke. Pinsonneault fc Sills] 
||2004) . However, a detailed comparison is difficult because 
they adopted input physics different from ours and they 
recalibrated the mixing length parameter and initial 
helium abundance on the Sun for each perturbed set of 
models. A variation of the initial Y and/or qml in the 
perturbed models due to the solar re-calibration might 
partially counterbalance or increase the effect induced 
by the sole variation of the analysed quantity. For this 
reason, we preferred to show the contribution on the LDB 
age of the sole perturbed quantity with all the other 
parameters fixed. Moreover, the use of a solar calibrated 
Qml does not guarantee a better agreement with radii of 
low pre-MS stars, where the lithium is actually depleted. 
Indeed, several papers have shown that in these stars the 
superadiabatic convection is much less efficient than in the 
Sun and that a proper value of qml is of the order of 1 
(see e.g. Venjmra^^t_^ I ll998l . iD’Antona fc Montalban|l2003l . 
Landi n et al.l |2006|. iGennaro . Prada Moron i fc Tognellil 


2012I . iTognelli. Degi lnnocenti fc Prada Moronil 12012 . and 
references therein). 


3.1 Nuclear cross-sections 

Since we are interested in LDB, the only nuclear reactions 
that might have an effect are those that take place before 


Table 3. Input physics varied in the computation of perturbed 
stellar models and their assumed uncertainty or range of variation 
(see text). The flag ‘yes’ specifies the quantities taken into account 
in the cumulative uncertainty calculation (see Sect. [Hi. 


Quantity 

Error 

Global 

^H(p,7)^He reaction rate 

±3% 

No 

reaction rate 

±5% 

No 

•^H('^H,p)'^H reaction rate 

±5% 

No 

^Li(p,Q!)ci reaction rate 

±10% 

Yes 

Electron screening(pH-^Li) 

±50%, ±100% 

No 

^rad 

±5% 

No 

Tph 

2/3, 100 

Yes 

BCs (“1 

AHFll, KS66 

No 

EQS 

OPAL06, FreeEOSOS, 
SCVH95 

No 


AHFll:lAllard Homeier fc Frevtad <1201111 : KS66: 

[Krishna Swamvl (|l966l) ' 

W OPAL06: iRog eiS-^Jjayfonoyl 1200^; FreeEO SOS: IrwinI 
ll2008h : SCVH95: ISaumoinChabrier^ van Hoi^ lll995 1 
The BCs/EOS have been varied by using tables provided by 
different authors, because a proper evaluation of the uncertainty 
is lacking. 



log l/l,„ 


Figure 2. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models and the sets with per¬ 
turbed ^Li(p,o)o cross-section. 


or during the ^Li-burning. Such reactions are: ^H(p, 7 )^He, 
^H(^H,p)®H, ^H(^H,n)®He, and ^Li(p,a)Q. 

We adopted the following uncert ainties in the quoted 
react ions: ±3% for the p+^H reaction dPescouvemont et al.l 
|2004| . for temperatures of about 1 0^ K, typical of d-bur ning), 
±5% for the ^H-|-^H channels_j[Tuminoet_a^ 2014l l. and 
±10% for the ^Li±p reaction dLamia et al.l 2012 1. 

Notwithstanding the crucial role played by deuterium 
burning during the early pre-MS evolution, due to the slight 
variation of the reaction rates within the current uncertain¬ 
ties, the effect on the LDB age is completely negligible. For 
this reason we do not show the related plots. Fig. [5] shows 
the relative age difference, i.e. (perturbed model — reference 
model) /reference model versus luminosity, due to the uncer¬ 
tainty in the ^Li(p,a)Q cross-section. The relative age differ¬ 
ence, about ±1%, is independent of the mass over the whole 
selected mass range. An increase of the cross section results 
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Figure 3. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models and the sets with the 
labelled variations of the electron screening factor. 


in a higher ^Li-burning efficiency at a given temperature, 
which leads to a more rapid ^Li depletion and to a lower 
LDB age at a fixed luminosity. This effect, analysed here for 
the first time, is small but systematic. We computed also ref¬ 
erence and perturbed models adopting the solar-calibrated 
mixing length parameter (i.e. aML=1.74), verifying that the 
relative LDB age differences are completely unaffected by a 
variation of aML- 


3.2 Plasma electron screening 


An important point to discuss is the effect of the elec¬ 
tron screening on ^Li-burning. Plasma electrons around 
the interacting nuclei reduce the effective Coulombian 
repulsion enhanc i ng the reaction rate by a factor /d 
(see e.g. Salget^ IfiSdI . iDewitt, Graboske fc Cood^ Il97d . 
iGraboske et al. 1973ll . A similar effect, due to atomic elec¬ 
trons, is present in the measurements performed in the 
laboratory. To this regard, there are hints that atomic 
electron screening measured in laboratory is systematically 
lower (^ 1/21 tha n the theoretical expectations (see e.g. 
IPizzone et al.llioid . and references therein). 

It is not yet clear whether theoretical computations of 
plasma e lectron screening are a ffected by a similar problem 
(see e.g. ICastellani et al.l ll996l ). Moreover, at the moment 
a reliable estimate of such an uncertainty source is lacking. 
However, the effect of modifying the plasma screening dir¬ 
ectly affects the reaction rate. For what concerns the LDB 
age, the uncertainty in plasma screening introduces a sys¬ 
tematic effect on the central temperature value at which 
a reaction becomes efficient, thus reflecting eventually on 
the LDB age. Given such a situation, it might be useful to 
check the effect of a variation of /pi in the ^Li(p,a)Q reac¬ 
tion. Given the lack of a solid estimate of the uncertainty in 
/pi(Li), the best thing we can do is assuming a maximum 
discrepancy between the predicted and real electron plasma 
screening efficiency similar to that observed in the laborat¬ 
ory; thus, we computed two sets of perturbed models with 
/pi(Li) increased by a factor 1.5 and 2.0. 

Fig. 0 shows the relative age difference versus lumin¬ 
osity due to the variation of the ^Li(p,a)Q plasma electron 
screening. Increasing the electron screening leads to higher 


reaction rates and hence to lower LDB ages at a given lumin¬ 
osity. An enhancement by a factor of 1.5 and 2.0 produces, 
respectively, a relative age decrease of the order of 3% and 
5.5% at a fixed luminosity. These results do not change if 
Qiml=1.74 is adopted in place of aML=1.00 to compute both 
the reference and perturbed sets of models. 

The plasma electron screening variation, although never 
analysed before, has a significant effect on the LDB age. How¬ 
ever, there are not yet robust estimates of the uncertainty 
in the ^Li(p,a)Q: plasma screening factor, so present com¬ 
putations are intended to only give an idea of this effect 
on the LDB. For this reason, we decided to not include this 
quantity in the cumulative uncertainty computed in Sect. 
[5] with the warning that a possible uncertainty source has 
been neglected. 


3.3 Outer BCs 


To solve the differential equations which describe stellar in¬ 
teriors, suitable outer BCs are required. The common ap¬ 
proach consists in specifying pressure P(rph) and temperat¬ 
ure r(rph) provided by a detailed atmospheric model at a 
given optical depth Tph (the matching point between atmo¬ 
sphere and interior). To estimate the effect of BCs on LDB 
age, two aspects must be considered: (a) the adopted atmo¬ 
spheric model and (b) the choice of Tph. 

The adoption of a detailed non-grey atmospheric struc¬ 
ture becomes important when the star has a thick con¬ 
vective envelope and the temperature profile is sensitive to 


the heat tran s port in the thin atmospheric layers (see e.g. 
Aamaj][j]T69|^jDOTmmu Nelgon_j^^]hm^^ et al.l 


1997I. Baraffe et al.ll 19981. Montalban et al.1 20041. 

Unfortunately, detailed atmospheric model tables do 
not contain the related uncertainties. Given such a situ¬ 
ation, one possibility to address the impact of the adopted 
BCs on LDB age is to compute models with BCs obtained 
from different atmospheric calculations. In the following we 
analysed the effect of the adopti on of two detailed non-grey 


anaiysea tne enect ot tne aaopti on ot two detaiiea non-grey 
atmospheric tables, namely the Brot^^^^auschild^_|200^ 
BH05) , our reference, and the lAllard. Homeier fc Frevtad 
(l201ll . AHFll), computed by me ans of the same hydrostati c 
atmospheric code (PHOENIX: [Hauschildt fc Baronlll999l L 
In both cases we used Tph = 10 - We c omputed also models 
with the grey [Krishna Swamvl lll966l . KS66) T — r profile, 
and Tph = 2/3. We emphasize that the last set should be 
considered as an extreme case, since it is well known that 
the adoption of grey atmospheric model is only a rough ap¬ 
proximation of the atmo spheric structure of extremely cold 
low-mass stars (see e.g. iBaraffe et al.lfl9^ . and references 
therein). 

A further complication is given by the fact that the 
AHFll tables for low-mass stars are available for only one 
metallicity (Z = 0.0134). We decided to perform the com¬ 
parison at a fixed Z, so even the BH05 (our reference) and 
the KS66 models have been computed for such a metallicity, 
which, only in this case, is different from the reference one. 

Fig.[3]shows the relative age difference versus luminosity 
due to the adoption of the quoted atmospheric models, for 
two ctml values. As expected given the similarity of the two 
atmospheric models, the effect of using AHFll in place of 
BH05 is relatively small, being of the order of 1% over a 
large interval of luminosity and reaching a maximum value 
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Figure 4. Relative age difference at LDB as a function of lumin¬ 
osity between models with the reference (BH05) and the labelled 
outer BCs, for aML= 1-00 (upper panel) and aML= 1-74 (bottom 
panel). As discussed in the text, in this case the metallicity value 
Z = 0.0134 has been used instead of the reference one. 

of about 3% only at very low luminosities. Note also that, in 
this case, the relative LDB age differences are essentially the 
same using q:ml=1.00 {upper panel) or the solar-calibrated 
value aML = 1.74 {bottom panel). 

A much stronger effect is due to the adoption of the 
grey KS66, which leads to relative age differences as large as 
8-15% for aML=1.00. These differences slightly reduce in the 
case of qml = 1-74 towards the upper bound of the luminosity 
for those masses where superadiabatic convection occurs in 
a progressively larger and larger region. 

Besides t he _ chosen atmospheric model, 

iMontalban et all (l2004l ) showed that the choice of Tph 
is important when atmospheric and interior computa¬ 
tions do not use the same input physics/prescriptions. 
In this sense, Tph is another free-parameter that usu- 
ally lies in the range [2/3, 1001 (see e.g. Table 2 in 
iTognclli. Prada Moroni fc Deg^Innocentill2011^ . 

Fig. [5] shows the relative age difference versus luminos¬ 
ity obtained varying Tph and keeping fixed the atmospheric 
model, i.e. BH05, for two mixing length parameter values. 
Models are computed for Tph = 2/3, 10 (reference), and 100. 
Passing from Tph = 10 to 100 the age at the LDB decreases 
whereas from Tph = 10 to 2/3 it increases; the largest differ¬ 
ence occurs when Tph = 2/3 is adopted. The extent of the 
relative age difference depends on the used mixing length 
parameter, the effect for oml^I-OO being roughly twice that 
for qml = 1-74. Note that all the atmospheric models ad¬ 



log i-/L,„ 


Figure 5. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models (Tph = 10) and the sets 
with the labelled Tph values, for Q!ml= 1-00 (upper panel) and 
aML= 1-74 (bottom panel). 

opt the same mixing length value (i.e. ctML.atm = 2.00, for 
r < "Tph). We are forced to keep it fixed because we do not 
have the code used to compute the atmospheric structure. 

In convective regions, the temperature gradient be¬ 
comes progressively more sensitive to the mixing length 
parameter as superadiabaticity gets larger. The models with 
Tph = 2/3 have a superadiabatic zone in the interior (t ^ 
Tph), where we can actually change oml, larger than the 
model with Tph = 100. The latter set of models has con¬ 
sequently a reduced sensitivity to the adopted mixing length 
parameter in the interior. 

Given the difficulty in estimating the uncertainty in the 
adopted BCs, for the computation of the cumulative error 
stripe in LDB age in Sect. [S] we took into account only the 
effect of changing Tph. 

3.4 Radiative opacity 

Although radiative opacity is one of the main ingredients in 
stellar computations, the current generation of tables do not 
contain any indication about the uncertainty. As a first step, 
to give an estimate of the opacity uncertainty propagation 
on the LDB age, we analysed the effect of an uncertainty 
in the Rosseland co efficients {K r ad) of ±5% in the whole 
structure, as done in I Valle et al.l ll2013al ). 

In convective regions, the actual temperature gradi¬ 
ent depends on the radiative opacity only in superadiabatic 
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Figure 6. Relative age difference at LDB as a function of lumin¬ 
osity between models computed with the reference and the ±5% 
perturbed radiative opacity. Grey outer BCs are adopted instead 
of the reference ones. 


Figure 7. Relative age difference at LDB as a function of luminos¬ 
ity between the reference models and those computed with radi¬ 
ative opacity tables extrapolated with different procedures. Grey 
outer boundary conditions are adopted instead of the reference 
ones. 


zones, while elsewhere it is mainly determined by the EOS. 
In the inner region of a fully convective star, given the high 
density, the temperature gradient is essentially the adia¬ 
batic one, while it becomes progressively more and more 
superadiabatic moving towards the surface, where the dens¬ 
ity drops. It is thus clear that in the interior of stars, the 
change of the opacity coefficients has a negligible effect on 
the structure, whereas the effect increases towards the sur¬ 
face. 

For models computed with detailed non-grey BCs, we 
can not modify the atmospheric opacity because we do not 
have the atmospheric code, thus preventing the possibility 
to consistently check the effect of Krad- Moreover, for the 
selected mass range in the interior (r ^ Tph), where we can 
actually modify the opacity coefficients, the density is high 
enough to guarantee an almost adiabatic convection. For this 
set of non-grey models we verified that the LDB age is not 
affected by a variation of the radiative opacity coefficients 
in the interiors (for both Tph = 2/3 and 10). 

A way in which we can estimate the effect of an opa¬ 
city variation extended also to the atmosphere is to com¬ 
pute stellar models adopting grey outer BCs. In this case 
we can actually calculate the atmospheric structure consist¬ 
ently with the Krad variation used in the interior. Although 
we are aware that the adoption of a grey BC is not the best 
choice for the mass range we are dealing with, such an ap¬ 
proach allowed us to give an estimate of the effect on LDB age 
of an opacity variation in the whole structure, atmosphere 
included. 

Fig. [6] shows the relative age difference versus luminos¬ 
ity caused by a variation of ±5% of the radiative opacity 
coefficients. Due to the progressively more extended super¬ 
adiabatic zone in the outermost layers at increasing stellar 
mass, the relative age uncertainty increases with the lumin¬ 
osity from negligible values at the faint end to about 1% at 
the bright end. At a given luminosity, the lower the opacity 
and the larger the central temperature and, consequently, 
the lower the LDB age. 

A further uncertainty source related to radiative opa¬ 
city is worth to be discussed. The radiative opacity tables do 
not cover the entire temperature—density plane required for 
computing pre-MS evolution for masses in the range [0.06, 


0.4] M fT) and an extrapolation to higher densities is nee ded. 
As in iTognelli, Prada Moroni fc Degl’Innocentil (1201 il l, in 
order to evaluate the impact of the extrapolation on the LDB 
age, we analysed three different extrapolation techniques: 
constant, linear from the last two points, and a linear fit 
from the last four points (our reference). As in the previous 
cases, we cannot properly evaluate the effect of changing 
the extrapolation on the LDB age for non-grey models, be¬ 
cause we can not modify the atmospheric models. The only 
thing that we can analyse for non-grey models is the effect 
of changing the extrapolation in the interiors. As expected 
the effect of such a change is negligible. 

However, to have an idea of the extrapolation impact 
we computed three additional sets of models with grey BCs 
where we can actually change the extrapolated opacities also 
in the atmosphere. Fig. [3 shows the relative LDB age differ¬ 
ence between models computed with the reference, i.e. four 
point linear, and the constant and two points linear radiative 
opacity extrapolations. The last two extrapolations lead to 
underestimate the LDB age with respect to the reference one; 
however, such an effect vanishes at high luminosities. The 
constant extrapolation produces an effect which is about 
three times larger than the linear one. 

A final uncertainty source related to radiative opa¬ 
city is the heavy elements distribution adopted to com¬ 
pute Rosseland coefficients. To analyse the effect on the 
LDB age of changing the heavy element mixture at a fixed 
total metallicity we computed models adopting the 
ICrevesse fc Sauv3 lll998l . GS98) solar mixture. As in the 
previous cases, we discuss only the effect on grey models 
since non-grey ones are completely unaffected also by the 
mixture change in the opacity tables used to compute the 
interiors, the mixture adopted by BH05 atmospheric mod¬ 
els being hxed. Fig. [8] shows the relative LDB age differ¬ 
ence between models adopting the reference (AS09) and the 
GS98 solar mixture and computed with the grey BCs. De¬ 
pending on the mass, the stellar structure crosses a region 
in the density—temperature plane where the opacity coef- 
hcie nts are particularly sensitive to different elements (see 
e.g. ISestito et al.|[200a l. Thus, the opacity variation due to 
the adopted mixture is a complex function of temperature 
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Figure 8. Relative age difference at LDB as a function of lumin¬ 
osity between models computed with radiative opacity adopting 
AS09 (reference) and GS98 solar mixtures. Grey outer BCs are ad¬ 
opted instead of the reference ones. 
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Figure 9. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models (OPAL06+SCVH95) 
and the sets with the labelled EOS. 


and density. However, as clearly visible in Fig. |8l the effect 
induced by the adopted mixture is very small, being always 
in the range [—0.5%, -1-0.5%]. 

We computed sets of models with the reference and per¬ 
turbed radiative opacity coefficients (Tcrad i 5%, different 
high-density extrapolations, different heavy elements mix¬ 
tures) also for q:ml = 1.74. The effect of an opacity vari¬ 
ation on LDB age is completely insensitive to a mixing length 
change. 

Summarizing the results of this section, we showed that 
changing the radiative opacity only in the interiors does not 
affect the LDB age. On the contrary, a sizeable effect is ex¬ 
pected when the radiative opacity variation is extended to 
the atmosphere. However, since we were able to compute 
such an effect only in the case of grey BCs, we preferred to 
neglect the opacity contribution in the computation of the 
cumulative uncertainty in the LDB age (see Sect. 0, being 
aware that a possible uncertainty source is missing in the 
cumulative error estimate. 


3.5 Equation of state 

The structure of fully convective pre- M S stars depends 


ine structure ot tully convective pre- M s stars depends 
on the EDS (see e.g., Mazziteii|_ 198^ D’Anton3 1 19931 . 


iTognelli. Prada Morimi fc Degrinnocentii|201lf). 

As for radiative opacity, the most commonly adop¬ 
ted EOS tables do not contain the uncertainty associated 
with the thermodynamic quantities (i.e. density, specific 
heat, adiabatic gradient etc.). Moreover, these quantities 
are strictly correlated among each other, thus preventing 
to follow a similar approach to estimate the uncertainty 
propagation into LDB age as that used for opacity. For 
this reason, to roughly estimate the dependence of the LDB 
age on the EOS, in addition to our reference set of models 
(OPAL06-I-SCVH95 EOS), we computed three additional 
sets adopting differe nt EOS table s. We used the sole OPAL06 
EOS, the FreeEOS (Ilrwinll2008l . in the EOSl configuration), 
and the sole SCVH95 EOS. Since both the FreeEOS and 
OPAL06 EOS do not cover the whole temperature—density 
plane suitable for computing the entire mass range [0.06, 
0.4] M© studied in this paper, we limit the comparison to 
the common mass interval. 


Fig. [9] shows the relative age difference versus lu¬ 
minosity due to the adoption of the quoted EOS. Both 
FreeEOS and OPAL06 induce age differences less than 
±1%. A larger effect 1-3% i s caused by the SCVH95 EOS. 
iBurke, Pinsonneault fc Silld (l2004ll already studied the ef¬ 
fect on LDB of the EOS, but they used tables different from 
those adopted here, thus preventing a detailed compar¬ 
ison. However, their results are, at least qualitatively, in 
agreement with ours. We also verified that the adoption of 
Qiml=1.74 in the reference and perturbed EDS models leads 
to the same relative LDB age differences. 

Due to the difficulty to properly estimate the uncer¬ 
tainty propagation of the adopted EOS, we did not take such 
an effect into account for the computation of the cumulative 
error bars in LDB age in Sect. [5] 


3.6 Mixing length efficiency 


Although commonly adopted, the mixing length theory 
is not a fully consistent treatment of the superadiabatic 
convection and, consequently, it introduces a relevant 
uncertainty source in computation of stellar structures with 
extended convective envelopes. The commonly adopted 
solar calibration of the mixing length parameter oml does 
not in principle guarantee a correct convective efficiency 
for stars, as in our case, of mass and/or evolutionary 
phase different from that of the Sun. In order to check 
the effect of changing the mixing length value on the 
LDB age, we computed two sets of models with aML = 
1.74 (our solar calibrated) and the other with oml^ 
1.00, as sugg e sted by some authors for pre-MS stars 
(see IStassunf l2008l. iGennaro, Prada Moroni fc Tognellil 
2 OI 2 I. ITognelli. Degl’Innocenti fc Prada Moronil 20121 . 
Somers fc Pinsonneaullllioid 1 

Fig- US] shows the relative age difference versus lumin¬ 
osity due to the adoption of the quoted mixing length para¬ 
meters. The effect of changing oml from 1.00 to 1.74 is to 
increase the central temperature, hence to increase the effi¬ 
ciency of ^Li-burning, and to reduce the age at the LDB. The 
relative age difference of the LDB is of the order of about 
— 1% and —6% for the lowest and, respectively, highest lu¬ 
minosities. The effect is lower at very low luminosity (very 
low mass stars) as a result of the thinner superadiabatic 
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Figure 10. Relative age difference at LDB as a function of lu¬ 
minosity between sets computed with two different mixing length 
parameter values, oml= 1-00 (reference) and «ml= 1-74 (solar 
calibrated). 


region. Such a result is in a greem ent with that shown by 
iBtirke. Pinsonneault fc Sills! (120041 ). 


4 INDIVIDUAL CHEMICAL COMPOSITION 

UNCERTAINTIES 

Besides the input physics, pre-MS stellar models depend 
also on the adopted chemical composition, namely the ini¬ 
tial metallicity Z, helium Y and deuterium Xd abundances. 
Thus, when comparing observations and theoretical mod¬ 
els, the proper chemical composition should be adopted. 
However, quite often only independent measurements of the 
current photospheric [Fe/H] of the stellar cluster/group are 
available, whereas direct estimates of Y, Xd, and metal 
abundances are missing. To overcome such a lack various 
assumptions are needed. 

The [Fe/H] value is converted into the total metallicity 
Z adopting, for Population I stars, a solar-scaled distribu¬ 
tion of metals, and assuming an initial helium abundance. 
Stellar modellers usually compute initial Y using the follow¬ 
ing linear relation: 

Ay 

y = yp + (1) 

where Up is the primordial helium abundance and AYjAZ 
is the helium-to-metal enrichment ratio. The proper initial 
metallicity Z to be adopted in stellar models follows from: 

7 = C-Yp){Z/X)q 

10 -[Fe/Hl+ (1 +Ay/AZ)(Z/X) 0 ’ 

where {ZIX)q is the current solar photospheric metal-to- 
hydrogen ratio. For the refe rence set of m odels we adop- 
ted Yp = 0-2485 ± 0.0008 (ICvburtlliooi) , AF/AZ = 2 
llCasagrandel 1200^ . and {ZjX^Q — 0.0181 (lAsplund et al.l 
l2009l) . 

All the quantities in equations © and © are known 
with an error, shown in Table |31 which directly propagates 
into the final Y and Z values. Note that these errors are in¬ 
tended to represent the edges of the variability region, rather 
than confidence intervals. Similarly to what was done in the 
previous sections, as a first step we computed perturbed 
models by varying a single parameter (AF/AZ, [Fe/H], and 


Table 4. Chemical composition parameters varied in the compu¬ 
tation of perturbed stellar models and their assumed uncertainty 
(see text). The flag ‘yes’ in the last column specifies the quantities 
taken into account in the cumulative uncertainty calculation (see 

Sect. [Sj- 


Quantity 

Error 

Global 

[Fe/H] 

iO.ldex 

Yes 

AF/AZ 

±l 

Yes 

(Z/A)0 

±15% 

Yes 

Ad 

±l X 10-5 

Yes 
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Figure 11. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models (F = 0.274) and the 
sets with the labelled initial helium abundance Y. 


(ZIX)q) at a time keeping all the others fixed to the ref¬ 
erence value. In the following, we did not take into account 
the uncertainty in Fp, being negligible. 

To the best of our knowledge, this is the first detailed 
analysis of the chemical composition uncertainty in the LDB 
age estimates. 


4.1 Initial helium abundance 


Following equation ©, the initial helium abund¬ 
ance adopted in stellar models strongly depends 
on the helium-to-metal enrichment ratio, which 

is poorly constrained by observat ions _ (see e.g. 

iGennaro. Prada Moroni fc Degl’Innocentil 201^ We 


adopted an error on AF/AZ of ±1 I Casagrandd l2007l ). 


and we kept fixed the value of Fp. In order to quantify 
the effect of the quoted uncertainty on the LDB age, we 
computed two sets of models with the reference Z = 0.0130 
and Ad = 2 x 10“®, and two different values of F, namely 
0.261 and 0.287, which correspond to AF/AZ = 1 and 3. 

Fig. llll shows the relative age difference versus luminos¬ 
ity due to the adoption of the quoted initial helium abund¬ 
ances. The relative age uncertainty is roughly ±2% over the 
whole mass range. The larger the helium abundance, the 
higher the central temperature, and the earlier the lithium 
depletion at a given luminosity. The results do not change 
if aML=1.74 is used instead of qml=1.00. 
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Figure 12. Relative age difference at LDB as a function of lumin¬ 
osity between the reference set of models {Z = 0.0130) and the 
sets with the labelled initial metallicity Z. 
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Figure 13. Relative age difference at LDB as a function of lumin¬ 
osity between models adopting the reference {(ZIX)q = 0.0181) 
and the labelled {ZIX)q values. 


4.2 Initial metallicity 

In order to analyse the effect of the metallicity on the LDB 
age, we computed two sets of models with the reference 
Y = 0.274 and Xd = 2 x 10~®, and two different metallicity 
values, Z = 0.0105 and Z — 0.0155, which roughly corres¬ 
pond to the typical observational error in [Fe/H], namely 
± 0.1 dex. 

Fig. [T5] shows the relative age difference versus lumin¬ 
osity due to the adoption of the quoted initial metallicities. 
The effect is almost negligible at faint end of our range and 
it increases up to about ±2% at higher luminosities. Such a 
behaviour is qualitatively similar to that shown in Fig. as 
a variation in metallicity mainly translates into a variation 
in the radiative opacity. This is essentially the consequence 
of the metallicity dependence of the outer BCs. We verified 
that varying the metallicity only in the interiors, keeping 
fixed the BCs, produces a very small (< 0.5%) and constant 
effect on the LDB. In this case at a fixed luminosity, a larger 
Z leads to a larger mass (higher central temperatures) and 
consequently to a lower LDB age. Moreover, we verified that 
decreasing the luminosity, the pressure and temperature at 
the base of the atmosphere of models at the LDB get pro¬ 
gressively less and less sensitive to the metallicity. At higher 
luminosity, the effect of Z on the BCs becomes dominant 
and an increase of Z results in a larger LDB age. The com¬ 
parison between reference and perturbed models computed 
with q:ml = 1.74 gives the same results. 

4.3 Heavy elements mixture 

The value of (ZjX^Q in equation © depends on the ad¬ 
opted solar heavy elem ents mixture. I n our reference set 
of models we used the lAsplund et al.l (l2009l ) one, which 
leads to (Z/X)© = 0.0181. In order to quantify the im¬ 
pact of varying this quantity on LDB age, we computed two 
sets of models with the reference values of XYjXZ = 2 , 
[Fe/H] = ±0.0, and Xd = 2 x 10~^, adopting an uncertainty 
of ± 15% in (Z/X)( 7 ) Bahcalh Serenelli fc Pinsonneaulu 
I 2 OO 4 I . iBahcall fc Serenelli I 2 OO 5 I ). This is also roughly rep- 
resentative of the range of values spanned by differ¬ 
ent, but still l argely used, solar mixtures ((ZIX)i^ ~ 
0.0165 - 0.0244: iGrevesse fc Noelslll993l . [&evesse &: Sauvall 


1 19981 . lAsplund. Grevesse fc Sauvai2005l . ICaffau et al.ll200^ ) . 

Hence, using the lower and upper values of (Z/X)q we 
computed the corresponding models with {Y, Z) = (0.271, 
0.0120) and (0.278, 0.0150). 

Fig-Hnishows the relative age difference versus luminos¬ 
ity due to the adoption of the quoted initial metallicities and 
helium abundances. The shown behaviour is the consequence 
of a variation of both Z and Y resulting from a change of 
{Z/X)q at a fixed [Fe/H] and XYjXZ. Note that, as shown 
in Figs. [ 11 ] and 1121 increasing the luminosity, the effect of 
a variation on Z is opposite and partially counterbalanced 
by that of varying Y, leading to a total effect lower than 
±1%. The relative LDB age difference between reference and 
perturbed models is not affected by changing oml from 1.00 
to 1.74. 

4.4 Initial deuterium abundance 

Although its tiny initial abundance, deuterium plays a rel¬ 
evant role in pre-MS evolution because the energy released 
during d-burning temporarily slows down the gravitational 
contraction. In order to evaluate the effect of the initial deu¬ 
terium abundance uncertainty on the LDB age, we computed 
two sets of models with Xd = 1 x 10“® and 3 x 10“® repres¬ 
entative of the current Xd range in t he Galactic disc /local 
bubble/solar neighbourhood (see e.g.. ISembach|[201^ . 

Fig. [14] shows the relative age difference versus lumin¬ 
osity due to the adoption of the quoted initial deuterium 
abundances. The effect is always less than ±1%. The larger 
is the Xd, the longer is the d-burning phase, and thus the 
later is the lithium depletion at a given luminosity. The rel¬ 
ative LDB age difference due to the initial deuterium abund¬ 
ance variation does not change if q;ml=1.74 is used. 


5 CUMULATIVE PHYSICAL AND CHEMICAL 
UNCERTAINTIES 

In the previous section, we analysed the effect on LDB age 
of varying a single input physics at a time by keeping fixed 
all the others, and varying the chemical composition for the 
refe rence input physics. However, although commonly adop - 

ted llBildsten et aI.l[l997l . lBurke. Pinsomieault fc Sillal2004 ). 





























Cumulative uncertainties in LDB age 11 



log L/L 


Figure 14. Relative age difference at LDB as a function of lumin¬ 
osity between models adopting the reference {X^ = 2 X 10“®) and 
the labelled initial deuterium abundances. 

such an approach does not take into acconnt possible inter¬ 
actions among the different ingredients of stellar models. In 
other words, it assumes that non-linear effects due to simul¬ 
taneous variation of snch quantities are absent or negligible. 

For this reason, we decided to fol low a more snitabl e 
procednre similar to that described in IValle et al.l (l2013al ) , 
consisting in the systematic and simultaneous variation of 
the main input on a fixed grid. More in detail, stellar mod¬ 
els and hence the resulting LDB age depend on a vector of 
parameters {pj} (i.e. input physics and chemical elements 
abundances, as discussed in the previous sections). Each 
parameter pj can assume three values, namely pj,o — Spj, 
Pjfi, and pjfi + 5pj, where pjfl represents the reference value 
and 5pj the adopted uncertainty. In order to cover the whole 
parameter space of the simultaneously varied input, we com¬ 
puted a set of perturbed stellar models for each possible 
combination of the vector {pj}. Such a distribution-free tech¬ 
nique relies only on the specification of a sensible range of 
variation for each input and it does not require an explicit 
specification of the parent distributions of the varied para¬ 
meters. Thus, it allows us to determine the edges of the 
variability region rather than a confidence interval. Such a 
technique is particularly useful when no information about 
the errors distribution on the analysed quantity is available. 

This method is more robust than assuming an a priori 
independence of the input physics/chemical parameters im¬ 
pact, but it has the serious disadvantage of being much more 
computationally time consuming than the classic approach. 
So, in order to avoid a waste of time in useless computa¬ 
tions, we did not vary all the input physics studied in Sect. 
[3] We did not perturb the d-burning cross-sections (p+^H, 
since we showed that they have a completely neg¬ 
ligible effect on the LDB age. Moreover, as discussed in the 
previous section, we did not vary the EOS, radiative opacity 
and atmospheric models since we can not properly take into 
account their effects. The list of the input physics allowed to 
vary is shown in the last column of Tables [3] and [I] Note that 
the perturbed stellar models can not be ruled out by the cur¬ 
rent observational constraints. In fact, we checked that the 
relative radius variation due to the perturbed input physics 
at fixed mass in the very low mass regime is smaller than 
1 %. 

Regarding the chemical composition, the number of {Y, 


Table 5. Pairs of {Y,Z) values adopted for the computations 
of the models with the perturbed chemical composition used for 
constructing the cumulative error stripe (see text). 



[Fe/H] 

AYjAZ 

{ZIX)q 

Y = 0.2740 

Z = 0.0130 

-1-0.0 

2 

AS09 

Y = 0.2790 

Z = 0.0100 

-0.1 

3 

AS09 

Y = 0.2650 

Z = 0.0160 

-bO.l 

1 

AS09 

Y = 0.2750 

Z = 0.0090 

-0.1 

3 

AS09-15% 

Y = 0.2670 

Z = 0.0190 

-bO.l 

1 

AS09-|-15% 

Y = 0.2570 

Z = 0.0090 

-0.1 

1 

AS09-15% 

Y = 0.3020 

Z = 0.0190 

-bO.l 

3 

AS09-|-15% 

Table 6. Pairs of {Y,Z) values adopted for computing models 
with the individual perturbation of [Fe/H], AYjAZ^ and {Z/X)q 
(see text). 

n. z) 

[Fe/H] 

AYjAZ 

{ZIX)q 

Y = 0.2710 

Z = 0.0120 

-1-0.0 

2 

AS09-15% 

Y = 0.2780 

Z = 0.0150 

-1-0.0 

2 

AS09-|-15% 

Y = 0.2620 

Z = 0.0130 

-1-0.0 

1 

AS09 

Y = 0.2870 

Z = 0.0130 

-1-0.0 

3 

AS09 

Y = 0.2690 

Z = 0.0100 

-0.1 

2 

AS09 

Y = 0.2810 

Z = 0.0160 

-bO.l 

2 

AS09 


Z) pairs resulting from all the possible combinations of 
[Fe/H]±0.1 dex, XYfXZ ± 1, and (^/X)© ± 15% values 
is 27 (= 3®). However, it is not needed to calculate set of 
models for each (Y, Z) pair. Indeed, being interested only in 
the edges of the variability region, it is enough to compute 
stellar models for the 7 (T, Z) pairs listed in Table O which 
correspond to the reference set and to the sets with the ex¬ 
treme combinations of helium abundance and/or metallicity. 

In summary, we have adopted for the cumulative er¬ 
ror evaluation seven pairs of (Y, Z), three initial deuterium 
abundances Xd, two mixing length values oml, three Tph 
and three ^Li(p,a)Q: cross-section values, for a total of 378 
(= 7x3x2x3x3) sets of models. For each set we computed 
a grid of 31 stellar masses, for a total of 11718 stellar tracks. 

From the perturbed stellar models, we then obtained 
the cumulative error stripe, following the algorithm of 
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Figure 15. Cumulative error stripe on the relative LDB age com¬ 
puted by taking into account the stellar models with simultan¬ 
eously perturbed input physics and chemical composition, for 
“ML= 1-00 (upper panel) and oml= 1-74 (bottom panel). Er¬ 
rors computed by linearly and quadratically adding uncertainties 
due to the independent variation of each input physics/chemical 
composition are overplotted (thick dashed and thin dot-dashed 
line, respectively). 

constructing the c onvex hull described in appendix A in 
IValle et al.l (l2013al ). Here we briefly summarize such a pro¬ 
cedure. Each set of perturbed stellar models provides an 
LDB curve in the age—luminosity plane. Thus, at a given 
luminosity an ensemble of perturbed LDB ages is available. 
The convex hull of such an ensemble provides the error in 
the predicted LDB age for the chosen luminositjQ. The error 
stripe is then obtained by changing the luminosity. 

Our main aim is to quantify the cumulative uncer¬ 
tainty resulting from input physics and chemical compos¬ 
ition, keeping fixed oml- However, the oml value is itself 
very uncertain and it might significantly affect the estim¬ 
ate of the cumulative uncertainty. For this reason. Fig. [15] 
shows the error stripe on the relative LDB age computed by 
taking into account the stellar models with simultaneously 
perturbed input physics and chemical composition for two 
fixed oml values, namely q:ml= 1-00 and 1.74. 

The cumulative error stripe is asymmetric and it gets 
progressively broader an broader at increasing luminosity. 
The case with oml^I-OO shows the largest and most asym- 

^ In this case the convex hull coincides with the maximum and 
minimum values of the LDB age at the chosen luminosity. 


metric stripe. The minimum relative LDB age uncertainty is 
of about ±5% while the maximum is of -1-15%/—8%. The 
case with q;ml=1-74 (solar calibrated) shows a thinner and 
more symmetric error stripe, with a maximum extension of 
about -1-10%/—8%. 

Both the stripe asymmetry and dependence on oml res¬ 
ult essentially from the variation of Tph, which assumes non 
symmetric values with respect to the reference one (see Sect. 
[Oil . Indeed, we verified that if the contribution of Tph is 
removed from the cumulative uncertainty estimate, the er¬ 
ror stripe becomes essentially symmetric and unaffected by 
oml- 

Fig-HSjalso shows the relative LDB age uncertainty com¬ 
puted by linearly/quadratically adding the errors due to the 
individual variation of each input physics/chemical compos¬ 
ition used to build the stripe. To do this we computed ad¬ 
ditional sets of models to account for the effect of the in¬ 
dividual variation of each parameter used to obtain Y and 
Z (i.e. [Fe/H], AYjAZ, and {Z/X)q), as listed in Table 
[6| Such models have been obtained by keeping fixed the 
input physics and varying a single chemical input within 
its uncertainty range (i.e. [Fe/H] ±0.1, AYjAZ ± 1, and 
[ZIX)q ± 15%) in equations (|T]) and (|21l. The comparison 
is intended to quantify the difference between the classical 
method (linear/quadratic individual error sum), often used 
in the literature, and the simultaneous variation on a fixed 
grid (cumulative error stripe) presented in this paper. It is 
evident that the adoption of a linear sum produces an un¬ 
certainty region much closer to the cumulative error stripe 
than that obtained by a quadratic sum. In particular the 
latter method gives a systematic underestimation of the un¬ 
certainty region (about two times smaller). Fig. [TJjhas been 
computed perturbing simultaneously the input physics and 
the chemical composition. However, to better understand 
the results it is worth to discuss the two effects separately. 

First, we verified that the contribution of the input 
physics (keeping the chemical composition fixed to the refer¬ 
ence va lue) sums linearly, as already obtained bv IValle et al.l 
ll2013ah . In other words, in the case of a constant chemical 
composition and small perturbations of the input physics, 
the uncertainty region obtained by a linear sum coincides 
with the edges of the cumulative error stripe, while the quad¬ 
ratic sum results in a significative underestimate. 

Secondly, to obtain the uncertainty region edges due 
to the perturbed chemical composition using both the lin¬ 
ear and the quadratic sum, we used the additional mod¬ 
els listed in Table jS] It is important to notice that equa¬ 
tions m and © are non-linear in the parameters [Fe/H], 
AY/AZ and {Z/X)q. In addition, even if a symmetric per¬ 
turbation of the parameters is adopted the resulting Y and 
Z variation is not necessarily symmetric. The non-linearity 
of eqs 0 and © introduces a non linear response of the 
LDB age variation that essentially depends on the initial 
helium abundance and metallicity perturbation {SY, SZ) 
with respect to the reference values. In the case of the 
simultaneous perturbation of all the chemical parameters 
the initial helium abundance and metallicity variations are 
simply given by 6Ycu. = Tref - Tcumuiative and SZcu. = 

^ref ^cumulative (F/umulative and .^cumulative are the ValueS 

listed in Table ©, while in the case of the individual para¬ 
meter perturbation by dTin. = [Fi-ef - Tindividuai((5j5j)] 
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Figure 16. As in Fig. IlSI but taking into account only chemical 
composition uncertainties. 


and [■^ref ■^individual(*^Pj)] (individual) and 

•^individual (<5pi) are the values listed in Table [ 6 ] and Spj is the 
jth-parameter variation). The non-linearity of equations ([Tjl 
and © leads to ((JITu., SZeu.) ^ (5Tin., 5Zin.), thus to a 
different variation of the LDB age if the cumulative or inde¬ 
pendent linear sum models are considered. We verified that 
this produces the differences between the linear sum edges 
and the cumulative stripe. 

Fig. [16] shows the cumulative error stripe obtained by 
varying the chemical composition parameters but keeping 
fixed the input physics to the reference ones, for oml = 1-00 
and 1.74. The uncertainty in the initial chemical composi¬ 
tion produces a not negligible variation of the LDB age, which 
is slightly asymmetric due to the asymmetric variation of Y 
and Z, regardless the adopted mixing length. The relative 
LDB age uncertainty ranges from about 3% to about 6 %, in¬ 
creasing at larger luminosities. When compared to the total 
error stripe (which includes also the input physics perturba¬ 
tion), the contribution of the sole uncertainty on the initial 
chemical composition is of the order of 40%. Fig. [12] also 
shows the linear/quadratic sum edges. As anticipated, the 
uncertainty obtained with a linear sum does not coincides 
with the error stripe even if the differences are quite small 
(about 1 %), while the quadratic sum drastically underes¬ 
timates the uncertainty (up to 1.5—2 times smaller). 

One should be aware that the global error stripe shown 
in Fig. [TS] is probably an underestimate of the actual un¬ 
certainty in the LDB age since, as explained in the previous 
sections, some potential contributions have been neglected, 


Figure 17. Age at the LDB as a function of luminosity for our 
reference set of models, with Z = 0.0130, Y = 0.274, = 

2 X 10“®, for «ml= 1-00 (upper panel) and oml= 1-74 (bottom 
panel). The cumulative error stripe is overplotted as a shaded 
area. 

given the difficulty to properly account for them. However, 
assuming a linear behaviour also for these contributions, we 
can give a rough estimate to the total LDB age uncertainty 
taking into account the additional effect of the outer BCs, 
the EOS, and the radiative opacity. As shown in Sect. [3] an 
uncertainty of 1 % is due to the adoption of a different non¬ 
grey BC (Fig. ]4]|, 1 — 2% to the EOS (Fig. [9]l, and 1% to 
an opacity variation (Fig. [ 6 ]|. Thus, adding these contribu¬ 
tions the uncertainty gets larger, ranging from a minimum of 
about ±9% to a maximum of 4-19%/—12%, for oml = 1.00, 
or to a maximum of 4-14%/—12%, for oml = 1.74. It would 
be worth to have an estimate of the uncertainty affecting 
the plasma electron screening factor, which, as shown in 
Sect. 13.21 might significantly affect the LDB age. Moreover, 
as already mentioned, we focused on standard pre-MS mod¬ 
els that do not take into account some potentially important 
physical mechanism. In particular it has been recently shown 
that magnetic fields an d star spots might signifi c antly affect 
the L DB age estimate ll Jackson fc Jeffrie |2014| . iMalo et al.l 
l2014h . The systematic bias due to neglecting these physical 
processes might be comparable or even larger than the cu¬ 
mulative error shown in Fig. 1151 

Finally, Fig. [17] shows our reference LDB curve in the 
(logL/L 0 , age) plane with overplotted the cumulative error 
stripe, for the aML= 1.00 (upper panel) and aML=1.74 set of 
models (bottom panel). 

As a final comment, our cumulative LDB age un- 
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certainty is larger (about 5 0%) than that provided by 
iBurke, Pinsonneault fc Sills! ll2004h . However, a detailed 
comparison of the two results is not straightforward be¬ 
cause of both the different uncertainty sources taken into 
account and the adopted estimate method. As an example 
of the former point, we accounted for the uncertainty 
in the adopted chemical composit ion, which is lacking in 
Burke, Pinsonneault fc Sills! (l2004l . Concerning the latter 
point, they adopted as best estimate of the total uncer¬ 
tainty a quadratic summation rather than the cumulative 
stripe described here. 


6 CONCLUSIONS 


In this work, we performed a systematic analysis of 
the main theoretical uncertainties affecting the age at 
the LDB. With respect to oth e r studies in the liter¬ 
ature (e.g. Bildsten et al.l 19971, Ushomirsky et ahl 


Ijeffries fc Navloil l200ll . iBurke. Pinsonneault &i Silld 12004) . 
which analysed the effect of changing an input physics at 
the time by keeping fixed all the others, the main novelty of 
our approach consists in taking into account the simultan¬ 
eous variation of the main input/parameters without an a 
priori assumption of independence. More in detail, we com¬ 
puted sets of pre-MS models for all the possible combina¬ 
tions of the perturbations of the ingredients allowed to vary. 
Such an approach is much more robust in presence of inter¬ 
actions between the varied quantities, but it clearly requires 
the computation of a huge number of models. For this work, 
we computed about 12 000 pre-MS evolutionary tracks. 

Besides the method, we improved the analysis of the 
LDB age uncertainty by studying physical error sources never 
discussed before, such as the plasma electron screening and 
the nuclear cross-sections (^Li-|-p ^H-fp). A fur¬ 

ther improvement with respect to previous studies is the 
detailed computation of the uncertainties propagation due 
to the adopted initial chemical elements abundances, i.e. the 
total metallicity, the helium and deuterium abundances, and 
heavy elements mixture. 

The results of the uncertainty computations are shown 
in Fig. [T7] where the LDB age—luminosity curve with the 
cumulative error stripe for both aML= 1-00 and 1.74 has 
been plotted. To the best of our knowledge, this is the first 
time that such an e rror stripe has been shown. _ 

As discussed in [Burke. Pinsonneault fc Sillsl (|2004| ). the 
LDB age estimate is more accurate at faint luminosity. Our 
detailed computations confirm that the error stripe gets pro¬ 
gressively broader at increasing luminosities. Moreover, we 
also showed that the cumulative error is asymmetric and it 
depends on the adopted oml value. More in detail, the set 
of models with q:ml = 1.00 shows a larger and more asym¬ 
metric error stripe, with positive relative age errors being 
about twice the negative ones, ranging from « 5 to « 15%. 
The uncertainty range reduces to a maximum of about 10% 
if the solar-calibrated mixing length value (ctml = 1.74) is 
adopted. A further contribution of the order of 4% should 
be added to the global error to take into account the un¬ 
certainty in the EOS, outer BC and radiative opacity, not 
explicitly accounted for in the error stripe computation. Fi¬ 
nally, an additional uncertainty source potentially sizeable 
and worth to be studied in more detail is the electron screen¬ 


ing factor. We proved that increasing it of 50% (100%) leads 
to a 3% (5.5%) variation in the LDB age. 

We showed that the cumulative error stripe computed 
by simultaneously perturbing the input physics is in excel¬ 
lent agreement with the simple linear sum of the LDB relative 
age differences obtained by individually perturbing the in¬ 
put physics themselves. Regarding the chemical composition 
uncertainties, the agreement is worse but always better than 
1%. On the contrary, the use of a quadratic sum (as generally 
done in the literature) results in a systematic underestimate 
of the total LDB age uncertainty, which is about 1.5/2 times 
smaller than that obtained with the cnmulative error stripe. 

Another result of this work consists in quantifying, for 
the first time, the effect of the initial chemical elements 
abundances uncertainty on the LDB age. We showed that the 
contribution of the sole chemical composition uncertainty is 
not negligible ranging from ± 3 to ± 6%, thus accounting 
for at least ~ 40% of the total error bndget; this part of the 
error stripe is almost symmetric and independent of ctml 
value. 
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